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^ ! Abstract 
OO ■ 

£N| ■ We study the temporal aspects of quantum tunneling as manifested in time-of-arrival 

experiments, in which the detected particle tunnels through a potential barrier. In particular, 
we present a general method for constructing temporal probabilities in tunneling systems 
that (i) defines 'classical' time observables for quantum systems and (ii) applies to relativistic 
particles interacting through quantum fields. We show that the relevant probabilities are 
defined in terms of specific correlations functions of the quantum field associated to tunneling 
particles. We construct a probability distribution with respect to the time of particle detection 

' that contains all information about the temporal aspects of the tunneling process. In specific 

cases, this probability distribution leads to the definition of a delay time that, for parity- 
symmetric potentials, reduces to the phase time of Bohm and Wigner. We apply our results 
to piecewise constant potentials, by deriving the appropriate junction conditions on the points 
of discontinuity. For the double square potential, in particular, we demonstrate the existence 
of (at least) two physically relevant time parameters, the delay time and a decay rate that 
' describes the escape of particles trapped in the inter-barrier region. 
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1 Introduction 



^ I The issue of the time that a quantum particle takes to tunnel through a potential barrier has been 
studied since the early days of quantum mechanics [HE]. The search for an answer to this question 
has led to several different candidates for the tunneling time — for reviews, see Ref.[3] — rather than 
to a single expression derived unambiguously from first principles. Many existing definitions imply 
that tunneling times saturate in the opaque-barrier limit, thus suggesting superluminal speeds for 
particles traversing sufficiently long barriers, a phenomenon known as the Hartmann effect p[] . 

In this paper, we identify tunneling time through its association with quantum temporal ob- 
servables, i.e, time variables whose value can be determined in specific experiments. A quantum 
temporal observable is a random variable, so its determination requires the construction of a prob- 
ability density rather than the specification of a single number. Hence, any temporal parameters 
that characterize the tunneling process must be identified in terms of a probability distribution 
associated to physical measurements on the tunneling particles. In particular, we propose a defi- 
nition of tunneling time in terms of the time instant that a detector, located far from the barrier, 
records tunneling particles, i.e., that a definition in terms of time-of-arrival measurements. In our 
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opinion, this is the most natural operational definition of tunneling time, because a remote detec- 
tor does not interfere with the tunneling process — unlike existing proposed definitions of tunneling 
time in terms of external actions on the barrier region [51 |6] . 

A time-of-arrival experiment typically involves a particle source and a particle detector sep- 
arated by distance L at rest with respect to another. Source and detector are equipped with a 
pair of synchronized clocks. The time of arrival is the difference t between the clock readings of 
detection and emission respectively. This 'classical' definition of the time of arrival can be ex- 
tended to quantum particles using the Quantum Temporal Probabilities (QTP) method [13]. The 
main difference is that for quantum particles, the time difference t is a random variable that may 
take different values in different runs of the experiment. Thus, a quantum description involves the 
construction of a probability density P(L, t) for the time of arrival. When a potential barrier is 
placed along the line between emitter and detector, all information about the temporal aspects 
of tunneling (the tunneling time and any other physically relevant parameter) is contained in the 
probability density P(L, t) [7J — for other approaches to tunneling time in terms of time-of-arrival 
measurements, see Refs. [8| 191 [TO]. 

The construction of time-of-arrival probabilities associated to tunneling systems is a non-trivial 
task. In fact, the definition of time-of-arrival probabilities has run into ambiguities analogous to 
the ones in the definition of tunneling time — for reviews, see Refs. [HJ [J2]. The difficulty originates 
from the special role of time in quantum mechanics, where the time t appearing in Schrodinger's 
equation is an external parameter and not as an ordinary observable, like position or momentum. 
As a result, the squared modulus of the time-evolved wave-function, \ip(x,t)\ 2 , is not a density 
with respect to t, and, hence, it cannot serve as a definition for the required probabilities. 

In this work, we employ a general method for the construction of quantum temporal proba- 
bilities associated to any experimental configuration, that was developed in Ref. [13J. The QTP 
method incorporates the detector degrees of freedom into the quantum description, so that the 
temporal probabilities are always defined with respect to specific experimental set-ups. The key 
idea is to distinguish between the roles of time as a parameter to Schrodinger's equation and as a 
label of the causal ordering of events [UJ [T5] . This important distinction leads to the definition of 
quantum temporal observables. In particular, we identify the time of a detection event as a coarse- 
grained quasi-classical variable [161 E] associated to macroscopic records of observation. The time 
variables in QTP correspond to macroscopic observable magnitudes, such as the coincidence of a 
detector 'click' with the reading of a clock external to the system. 

A key property of the QTP method is that it applies to any quantum system, including rela- 
tivists quantum fields. Besides time-of-arrival probabilities, the method has been applied to the 
modeling of particle detectors in high-energy processes [13] , and to the study of temporal correla- 
tions in accelerated detectors [18J. Earlier versions of QTP [19] have been employed in studies of 
non-relativistic tunneling times [7J, non-exponential decays [20J and time-extended measurements 

Using the QTP method, we derive a probability density P(L, t) describing time-of-arrival mea- 
surements on relativistic particles that tunnel through a potential barrier. The interactions of 
the particles with the barrier and the detector are described in terms of Quantum Field Theory 
(QFT). The definition of relativistic tunneling times in terms of QFT (rather than relativistic 
quantum mechanics [22]) is important, because a QFT treatment of tunneling times is essential 
for a resolution of the superluminality paradox [22] ■ The reason is that superluminal speeds imply 
a violation of the principle of local causality, and this principle is implemented in quantum theory 
only if the interactions are expressed in terms of local quantum fields [211 125] . 

The probability density P(L, t) derived here incorporates all temporal aspects of the tunneling 
process. For a specific class of potentials and initial states, the only timescale relevant to tunneling 
is a temporal variable id- This variable describes the delay in the transit time of particles crossing 
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through the barrier in comparison to the transit time of the same distance in absence of the barrier. 
For parity-symmetric potential, coincides with the Bohm-Wigner phase time [26J, obtained from 
the asymptotic analysis of wave-packets. However, our identification of the delay time does not 
employ non quantum-invariant notions, such as the peak of the wave packet or its center of mass, 
but follows from the structure of the probability density P(L,t), in specific regimes [27]. 

In general, a single time parameter, such as t d , does not capture all temporal aspects of the 
tunneling process. We demonstrate this point by studying tunneling through a double square bar- 
rier. We find that the probability distribution P(L, t) is characterized by two distinct timescales, 
a delay time td, as described earlier, but also a decay time T _1 that determines the escape rate 
of particles trapped between the two barriers. Previous work on this system, led by the a priori 
assumption of a single-time scale governing tunneling, had led to a very different physical descrip- 
tion of the tunneling process, according to which the tunneling time is largely independent of the 
inter-barrier distance [28| [29] . 

The main results of the paper are the following. 

1. We present a general methodology for defining the time-of-arrival probabilities P(L,t) in 
tunneling set-ups that is valid for any potential barrier and for any method of particle 
detection. We derive the probability distribution P(L,t) for relativistic particles interacting 
through quantum fields. We show that P(L, t) is obtained from a suitable two-point function 
of the quantum field smeared with a kernel that depends on the physics of the detection 
process. 

2. For the standard case of a potential barrier model by a background classical field, we derive 
an simple expression for the probability distribution P(L, t) as a positive linear functional of 
the initial state. The probability distribution P(L, t) can be explicitly constructed from the 
knowledge of the transmission and reflection amplitudes of the potential. 

3. We find that the total transmission probability coincides with the square modulus of the 
transmission amplitude associated to the potential only if the potential is parity symmetric. 
In general, they differ. This is because the QFT propagator through which the probability 
distribution P(L, t) is defined involves contributions from both left-moving and right-moving 
eigenstates of the singe-particle Hamiltonian. 

4. We construct the probability distribution P(L,t) explicitly for piecewise constant potentials, 
by deriving the appropriate junction conditions on the points of discontinuity of the potential. 
For the double square potential, in particular, we demonstrate explicitly that a single time 
parameter does not suffice to capture all information about the temporal aspects of tunneling. 

In a companion paper, we analyze in detail the superluminality paradox in tunneling, in light of 
the present results [23J. 

The structure of the paper is the following. In Sec. 2, we formulate the quantum tunneling 
of particles through a barrier in a quantum field theory language. In Sec. 3, we derive a general 
expression for the probability density P(L, t) of detection of particles through a potential barrier. 
In Sec. 4, we construct explicitly the probability density P(L,t) for piecewise-constant potentials. 
In Sec. 5, we summarize our results and discuss their further applications. 

2 QFT formulation of tunneling 

In this section, we formulate tunneling of particles through a potential barrier in a quantum field 
theory language. In particular, we consider a complex scalar field (f>(x), corresponding to particles 
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of mass m and charge e, interacting with a background static electric field. The electromagnetic 
potential is A^(x) = (A Q (x),0), where A (x) differs from zero only in a spatial region D. For 
simplicity, we restrict our considerations to one spatial dimension, and we choose the origin of the 
coordinate system so that D = [—d/2, d/2], where d is the length of the barrier. 
The classical Lagrangian density for the scalar field <p(x) is 

£(0, 0*) = |0 - iVct>\ 2 - \8 X <P\ 2 - m 2 |0| 2 , (1) 

where V(x) = eA (x). 

The associated Hamiltonian density H is 

U = M 2 + \d x (j)\ 2 + m 2 |0| 2 + iV(tx<P - vr*0*), (2) 

where n and n* are conjugate field variables to and 0* respectively. 

In order to quantize the field, we first introduce two copies H\ and Hi of the Hilbert space 
C 2 (R,dx) of square integrable functions on the real line. Hi and Hi are associated to a single 
particle and to a single antiparticle, respectively. The Hilbert space of the field is defined as 
J 7 (Hi) © J 7 (Hi). T(H) stands for the bosonic Fock space associated to the Hilbert space H, 

?{H) = C © H © (H ® H) s © (H ® H ® H) s © . . . , (3) 

where S denotes symmetrization. 

The creation and annihilation operators d(x) and d\x) for particles and b(x) and ^(x) for anti- 
particles are standardly defined on the Fock space J- (Hi) ^^(Hi). They satisfy the commutation 
relations 



[a(x), a* (a/)] = S(x — x') [b(x),lr(x')] = 5(x — x'), (4) 

with all other commutators vanishing. 

The creation and annihilation operators are also expressed in their smeared form 



d(f) = J dxd(x)f(x), d\f) = J dxd\x)f(x) (5) 
b(f) = [ dxb(x)f(x), b\f) = [ dxl\x)f(x), (6) 



where / is a square-integrable function on R. 

In order to represent the Hamiltonian Eq. (j2J) as an operator on the Fock space J 7 (Hi) ^)J-'(Hi) 
we introduce the field operators <f>(x),<fy(x) and their conjugate momenta 7r(x), tt(x). The field 
operators are defined in their smeared form as 

*(/) := \ (*(h 1/2 f) + b\h Q l/2 f)) rt(f) := ^ (a(h 1 ' 2 f) - a\h\' 2 f)) (7) 



4>\f) ■= ^ (Hho 1/2 f) + tf(h~ 1/2 f)) n\f) := -L (h(hl /2 f) - tf(hl /2 f)) . (8) 

where ho = \/—d^. + m 2 is the Hamiltonian operator for a free relativistic particle, defined on the 
Hilbert space Hi = Hi. 
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We construct the Hamiltonian operator H for the field by substituting the field operators above 
into the classical Hamiltonian density ([2]). After normal ordering, 

h = JdM (at WAl (,,,oa W + iH*)h 2 ( x , x -)b( x )) . (9) 

In Eq. (Q, h\(x, x') are matrix elements in the position basis of a Hamiltonian operator hi on the 
Hilbert space Hi of a single particle and h 2 (x,x') are matrix elements in the position basis of a 
Hamiltonian operator h 2 on the Hilbert space Hi of a single antiparticle. The operators hi t 2 are 
defined as 



hi = h + V, h 2 = h -V (10) 



where 



V := \ (C W /2 + hTvht") = V + [[V, hl'\ h?>\ (11) 

is the potential operator, that includes a QFT correction term. The sign difference between hi 
and h 2 indicates the fact that a potential barrier for particles is a potential well for anti-particles. 

In what follows, we will consider a positive- valued potential V(x), such that the corresponding 
single-particle Hamiltonian hi has purely continuous spectrum. The eigenstates of hi are charac- 
terized by energies E > m and they have a double degeneracy. We denote the pair of eigenstates 
corresponding to the same value E of energy as f k + and where k = \JE 2 — m 2 > 0. The 
eigenstates fk+{ x ) correspond to only positive momentum flux at x = oo and they are of the form 

MX) = V^\ x >( < 12 > 

In Eq. (II 2p . and Rk are the usual transmission and reflection coefficients for a right-moving 
plane wave of momentum k. The eigenstates are orthogonal to fk+ and they satisfy 

fk-( x ) = 7= 1 5 [-w k f k+ (x) + / fc+ (-a:)] , (13) 



where Wfc = -(T^R). + T k R*k) is essentially the overlap between f k+ and its parity-transform. 
For parity-symmetric potentials, i.e., for potentials such that V(x) = V(—x), the coefficient Wk 
vanishes and f k -(x) = f k+ (—x). 

The Hamiltonian h 2 , describing anti-particles, has both continuous (E > m) and discrete 
spectrum (E < m). We denote the continuous-spectrum eigenstates as g k+ and g k -. They are 
similar in structure to the generalized eigenstates f k ± of hi. Their form is not relevant to the 
purposes of this paper. 

We denote the discrete-spectrum eigenstates of h 2 as g n (x), where the label n runs into a finite 
set of values. The eigenstates g n (x) decay exponentially outside the barrier region. 

. v _ / A n e^ x x < -d/2 , . 

9n{X> \ B n e~^ x x>d/2 ' 1 > 

for some positive constants A n and B n , and 

7n = ^m 2 - El (15) 

For parity-symmetric potentials, B n = P(n)A n , where P{n) is the parity of the eigenstate g n (x). 
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Next, we define the energy-basis creation and annihilation operators, 

a k ± = J dxa(x)fl ± (x) b k ± = J dxb(x) g* k± (x) b n = J dxb{x)g* n {x) (16) 
d\ ± = J dxd\x)f k ±(x) b\ ± = J dxP(x)g k ±(x) P n = J dxP (x)g Tl (x). (17) 
Then, the Heisenberg-picture field operators are 

+ poo fit- / \ 1 

$(x, t) = I ^ {a k ,*hAx)e- iEkt + b{ a 9lA x ) eiEkt ) + E ^w}^^ (18) 

+ poo jr, i 

*) = E I ^ (ftik.^we-^-* + d\ji^y Ekt ) + E 7^^« (a;)e "^*- (19) 



Eqs. ( JTHl) and ( JT9l) are the basis for the QFT treatment of tunneling time in the following 
section. 

3 Probabilities for tunneling time 

In this section, we derive a general expression for the detection probability P(L,t) for relativistic 
particles tunneling through a potential barrier. First, we present a brief review of the QTP method, 
and then we derive an equation that relates the detection probability to the correlation functions 
of a quantum field theory. Finally, we specialize to the scalar field system described in Sec. 2, and 
we express the probability P(L, t) in terms of the transmission and reflection coefficient associated 
to the potential barrier. 

3.1 Background 

Our approach is based on the general methodology for constructing time-of-arrival probabilities 
associated to general detectors, that was developed in Ref. [33] (the Quantum Temporal Proba- 
bilities method). First, we give a brief review of the QTP method, setting up the background for 
its application to the tunneling-time problem in Sees. 3.2 and 3.3. 

The key result of Ref. [13] is the derivation of a general formula for probabilities associated 
to the time of an event (such as particle detection) in a general quantum system. The event time 
t is treated as a coarse-grained, quasi-classical parameter associated to macroscopic records of 
observation. 

Let % be the Hilbert associated to the physical system under consideration. We incorporate the 
measurement device into the quantum description, so % describes the degrees of freedom of both 
the microscopic particles and the macroscopic measurement apparatus. We assume that % splits 
into two subspaces: H = H + © H- . The subspace H + describes the accessible states of the system 
given that a specific event is realized; the subspace %_ is the complement of H + . For example, 
if the quantum event under consideration is a detection of a particle by a macroscopic apparatus, 
the subspace % + corresponds to all accessible states of the apparatus given that a detection event 
has been recorded. We denote the projection operator onto % + as P and the projector onto %_ 
as Q := 1 - P. 

Once the transition has taken place, it is possible to measure the values of a microscopic 
variable through its correlation to a pointer variable of the measurement apparatus. We denote by 
P\ projection operators (or, more generally, positive operators) corresponding to different values A 
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of a physical magnitude that can be measured only if the quantum event under consideration has 
occurred. For example, when considering transitions associated to particle detection, the projectors 
Pa rnay be correlated to properties of the microscopic particle, such as position, momentum and 
energy. The set of projectors Pa is exclusive (P\Py = 0, if A 7^ A'). It is also exhaustive given that 
the event under consideration has occurred; i.e., J2\P>~ = P- 

We also assume that the system is initially (t = 0) prepared at a state \ip ) e % +1 and that the 
time evolution is governed by the self-adjoint Hamiltonian operator H. 

In Ref. [13], we derived the probability amplitude A, [ti, £2]) that, given an initial state \ipo), 
a transition occurs at some instant in the time interval [£1,^2] and a recorded value A is obtained 
for some observable: 



\ij ;X,[t 1 ,t 2 ]) = - i e-^ T / dtC(X,t)\ij ). (20) 

Jti 

where the class operator C(X,t) is defined as 

C{X,t) = e iAt P x HS t , (21) 

and 

S t = lim (Qe- iflt/N Q) N (22) 

N— >oo 

is the restriction of the propagator into %_ 

If [P, H] = 0, i.e., if the Hamiltonian evolution preserves the subspaces 'H±, then \ip ; A, t) = 0. 
For a Hamiltonian of the form H = H + Hj, where [Hq, P] = 0, and Hj a perturbing interaction, 
to leading order in the perturbation 

C(X, t) = e lAot P x Hie- iflot . (23) 

The benefit of Eq. (123]) is that it does not involve the restricted propagator St, which is difficult 
to compute in practice. 

The amplitude Eq. f[2"U|) squared defines the probability p(X, [ti,t 2 ]) that at some time in the 
interval [t 1 , ^2] a detection with outcome A occurred 

P(A,M 2 ]) := ^o;\[ti,t 2 ]\^o;\[ti,t2\) = f 2 dt Pdt'Tr[C(\,t)po&(\,t)], (24) 

Jti Jti 

where p = \ip )(ip \. 

However, the expression P(A, [ti, £2]) does not correspond to a well-defined probability measure, 
because it fails to satisfy the Kolmogorov additivity condition for probability measures 

P(A, M 3 ]) = P(A, M 2 ]) + P(A, [t 2 ,t 3 ]). (25) 

Eq. (I25p does not hold for generic choices of ti,t 2 and £3. However, in a macroscopic system (or 
in a system with a macroscopic component) one expects that Eq. (125]) holds with a good degree of 
approximation, given a sufficient degree of coarse-graining [TBI E] • Thus, if the time of transition 
is associated to a macroscopic measurement record, there exists a coarse-graining time-scale a, 
such that Eq. (1251) holds, for \t 2 — 1\\ » a and \H — t 2 \ » o. Then, Eq. fT24"|) does define a 
probability measure when restricted to intervals of size larger than a. 

Eq. fl24l) simplifies if there is a separation of time scales between the macroscopic scale of 
observation, the coarse-graining time scale and the characteristic timescales of the microscopic 
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system. In this case, it can be shown [T3] that the probabilities Eq. ( |24|) are obtained from the 
probability density 



drTr 



C{X,t+ T -)p Q C\\,t- T -) 



(26) 



The probability density Eq. (|2~61) is of the form Tr[p Il(X,t)], where 

fl(A,t) = J dTC\X,t- T -)C{X,t + T -). 



The positive operator 



n T (iv) 



dt / dXIl T (X,t), 



(27) 



(2f 



corresponds to the alternative iV that no detection took place in the time interval [0, oo). The 
positive operator Il T (iV) together with the positive operators Eq. (j2"7|) define a Positive-Operator- 
Valued Measure (POVM). The POVM Eq. (l2~7]) determines the probability density that a transition 
took place at time t, and that the outcome A for the value of an observable has been recorded. 



3.2 Temporal probabilities from field correlation functions 

Next, we employ Eq. f)26p for constructing probabilities associated to time-of-arrival measurements 
in a quantum field theoretic set-up. The system under consideration consists of a quantum field 
that describes microscopic particles and a macroscopic particle detector. Eventually, we will 
identify the quantum field with the scalar field 4>(x), Eq. ffT8l) . associated to tunneling, but the 
results obtained in Sec. 3.2 apply for any quantum field theory. With small modifications, the 
method applies also to spinor and vector fields. The derivation presented here is an adaptation of 
the method described in Ref. [T3] . 

The Hilbert space for the total system is the tensor product T <g> H a . The Hilbert space % a 
describes the apparatus's degrees of freedom and T is the Hilbert space space associated to the 
field 0. 

The Hamiltonian of the total system is the sum Hq <g> 1 + 1 <g> H a + Hj. H is the Hamiltonian 
that describes the quantum field 0, H a describes the dynamics of the apparatus and Hj is an 
interaction term. 

The requirements of Lorentz covariance and unitarity in quantum field theory imply that the 
interaction Hamiltonian Hj must be a local functional of the fields 0(x) and <f>'(x). Therefore, we 
consider an interaction Hamiltonian of the form 



Hj = dx 



0(x) g> J\x) + f (x) ® J(x) , (29) 

where J and are current operators on the Hilbert H a of the detector. The ineraction Hamiltonian 
Eq. fl29l corresponds to particle detection by absorption. 

We also assume an initial state |^ ) °f the detector such that 

J(x)\^ ) = 0. (30) 

This condition guarantees that the detector is sensitive to particles rather than anti-particles. To 
see this, note that the term $(x) <S> J(x) in the Hamiltonian Eq. (129]) corresponds to the processes 
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of either a particle being created at the measurement device or an anti-particle being annihilated. 
The second process is negligible if the initial state contains only particles at energies E < m. The 
first process is not desirable in a particle detector. Eq. (|30|) guarantees that the only transitions 
in the detector are caused by the interaction with microscopic particles. 

Next, we specify the macroscopic variables associated to particle detection. These correspond 
to degrees of freedom of the macroscopic apparatus and they are expressed in terms of the positive 
operators 1 <g> III on J 7 (Hi © Hi) ® T-L a , labeled by the values L of a macroscopic observable 
correlated to a particle's position in the laboratory frame of reference. The operators IIl are 
defined on H a and they satisfy the completeness relation J dLYli = P, where P is the projector 

onto the subspace H + . We choose fl^ so that VTIx corresponds to an unsharp Gaussian sampling 
of position at X — L, 



n^p^jE f dXe-^\X,a)(X,al (31) 

The index a in Eq. ( IBTj) refers to the degrees of freedom of the apparatus other than the pointer 
variable, and S corresponds to the spatial resolution of the detector. 

We place no restriction on the apparatus' Hamiltonian H a , except for the requirement that it 
commutes with the operator P: [H a , P] = 0. This implies that the self-evolution of the apparatus 
degrees of freedom does not lead to spurious detection records. It follows that [1 ® P, H s ® 1 + 1 ® 
H a ] = 0; hence, the class operators C(X,t) follow from Eq. 

Substituting into Eq. f|27|) . we obtain 



P(L,t) = J dr J dxdx'(ifj \(f)' t (x',t~T/2)(f)(x,t + T/2)\tfj Q )R(L,T,x,x') (32) 
where <fi(x,t) and 4> da " er (x,t) are the Heisenberg-picture fields, is the initial field state and 

R(L,r,x,x') = (y \j(x')^ L e i6aT ^ L fi(x)\%). (33) 

is a kernel that involves the detailed physics of the detector. 

The macroscopic record of a detection must be correlated with the locus of the interaction, in 
the sense that yU^J^ (x)\^>q) ~ if \L — x\ is much larger than the detector's resolution 5. This 
implies that x and x' in Eq. ( 1331) lie within a distance 5 from L, otherwise the kernel R(L, r, x, x') 
vanishes. Equivalently, the kernel R vanishes unless X = (x + x')/2 lies within distance of order 8 
from L, and \z\ is of order S. If S is much smaller than the distance between source and detector, 
we can approximately set L = (x + x')/2, so that 

R(L, r, x, x') ~ S(L - ^-)g(x - x', r), (34) 

This form for R(L, r, x, x') is consistent with several different detector models that have been 
constructed in Ref. [13]. The kernel g(z,r) depends on the detailed physics of the detector. It 
vanishes for values of \z\ larger than the position resolution 5 of the detector and \t\ larger than 
the decoherence time of the detector, respectively. 
Eq. (132]) becomes 

P(L,t)= ! dr ( dz(i; \i\L-^,t-^)4>(L + ^,t + ^o)9(z,r). (35) 
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It is important to observe that no assumption about the form of the Hamiltonian H , or about 
the initial state |^o) has been made for the derivation of Eq. fl3"21) or Eq. (|3"oT) . Eq. fl3"21) reveals a 
relation between detection probability and field correlation functions that applies to any quantum 
field theory. The probability of detection is determined by the two-point function of the fields that 
interact with the detectors. 



3.3 The probability distribution for the time of detection 

Next, we specialize to the case of the fields (f>(x, t), $(x,t) of Eqs. (TT8T fT9l). associated to a 



tunneling set-up. We consider a single-particle initial state for the field \ip ) = J dxa) (x)ijjQ(x)\0) , 
where |0) is the field vacuum and ipo( x ) a single-particle wave function concentrated on positive 
values of momentum. 

Substituting Eqs. (TT8T [19]) into Eq. (135|) . we obtain 

P(L,t) = P (L,t) + Pt(t) + P 2 (L,t), 



(36) 



where 



P (L,t) 



dr / dzg(z, r) 



P 1 (L,t) 



J dxA(L + \ 

J dx'A(L-^,x';T-^)Mx' 
dzg(z, t) ^2 



2 



>b;t + -)ip (x) 



dr 



dr 



dk 



-ikz+iE k r 



-oo (27r)(2£, 



z z 
dzg{z, t)D(L - -, L + - ; -r) 



and 



A(x, x'; t) 
D(x,x';t) 



dk 



E / f— 

Jo V^Ek 



'x')e- lEnt . 



(37) 
(38) 
(39) 

(40) 
(41) 



The terms Pi(L,t) and P2(L,t) are independent of the initial state of the particle. They corre- 
spond to 'false alarms', that is spurious detection events — such terms are generic in the theory of 
relativistic quantum measurements [3T]. They act as noise that obscures the signal term P Q (L,t). 
The Pi term in is independent on the location of the detector. It is also independent of the barrier, 
as it persists for V(x) = 0. Its contribution to the total probability is negligible for detector sizes 
much larger than the de-Broglie wave-length of particles. The term P 2 originates from the anti- 
particle bound states and it depends on the the position L of the detector. However, the function 
D 2 (x,x' ,t), Eq. (l4"Tj) decays exponentially outside the barrier region, with a rate given by j n , Eq. 
(1TBT) . Therefore, the term P 2 (L,t) is strongly suppressed for L » max n {7 ? ^ 1 }. The largest values 
of 7 n are of the order of the length \Jd~jm « d. 

For a detector far from the barrier (L » d), the contributions from the false- alarm terms Pi 
and P 2 becomes negligible, and P(L,t) = P Q (L,t). 



The term Po(L,t) involves the function A(x,x';t) } Eq. ( T4T1) . For x > 



2 ' 



X < 



, we obtain 



A(x, x; t) 



dk 



+e 



ik{x+x') R*k + (^fc ~ 



ik(x—x') 



(T fc + w k - 



-R k + w k T k 



1 



+ e 



-ikjx+x') ~ w k 



w k R* k ){R k -w k T k 



+ e 



\Wk\ 

-ik(x-x') k 



w k FLj 
\w k \ 2 



wt 



-iEi.t 



(42) 
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The term in the second line of Eq. (|42|) vanishes when A(x,x';t) acts on wave functions with 
support on positive values of momentum. The second term in the first line of Eq. (|42p is strongly 
suppressed for x ~ L and t > 0. 

Thus, for sufficiently large L, the probability density P(L, t) becomes 

f°° rlh f°° rile' 

P(L,t)= - / - f(k, k')A k At,Mk)Mk')e i{k ' k W^-* 5 **, (43) 



where 

T k — WkR k 
5 



A k = (44) 



/(*, k') = dr dzg(r, z)e^^—^\ (45) 



ipo{k) is the Fourier transform of the initial state, and E k = Vk 2 + m 2 . 

For an initial state with support on positive momenta, P(L, t) is strongly suppressed for t < 0. 
Hence, we approximate J °° P(L, t) ~ f P(L, t). For a free particle (A k = 1), the time-integrated 
probability for a particle to be detected at x = L is 

'<M=jT (46) 

The quantity a(k) = then defines the absorption coefficient of the detector for particles with 
momentum k, i.e., the fraction of the number of incident particles absorbed by the detector. The 
function f(k, k') depends on k and k' only through the combinations (k + k')/2 and (E k + E k >)/2. 
For an initial state with momentum spread 8k much smaller than the mean momentum k, we 
approximate 

f(k, k') = ^\j{k~kyw\f(k',k')\. (47) 

Then, Eq. (l4"3l becomes 

P(L,t) = \A(L,t)\ 2 . (48) 



where 



dk 



A(L,t)= I ^a(k)\v k \A k ^ Q (k)e ikL - iEkt (49) 
Jo 2tt 

Eqs. pST — f4"9"|) are the main result of this paper. They express the probability distribution for 
particle detection in a tunneling experiments in terms of the transmission and reflection coefficients 
of the barrier, encoded in the coefficient A k . All information about the detector is encoded in the 
absorption coefficient a(k) for particles of momentum k and all information about the barrier in 
the complex amplitude A k . 

In absence of the potential barrier A k = 1, and Eq. fT45|) reduces to the time-of-arrival distri- 
bution for free relativistic particles derived in Ref. [13J, which generalizes Kijowki's probability 
distribution for the time of arrival of non- relativistic particles |30j . 

The time-integrated probability Eq. (HHj) is 

/dk 
a(k)\A k \ 2 \Mk)\ 2 . (50) 



11 



Eq. ( l5Uj) implies that \A k \ 2 is the transmission probability, i.e., the probability that a particle of 
momentum k crosses the barrier at any time t. For parity symmetric potentials, = T k and the 
transmission probability coincides with the standard expression \Tk\ 2 . In general, however, they 
differ. The reason is that in the QFT description of tunneling, both right-moving and left-moving 
modes contribute to field correlation functions, and by virtue of Eq. (135|) . they contribute to the 
detection probability. 

It is convenient to parameterize T k = \T k \e l< ^ k and R k = —i\Rk\e t ^ k+lXk . In parity-symmetric 
potentials, \k = 0. In the regime of opaque barrier, \T^\ « \R k \ — 1, hence, to leading order in 
\T k \ 

A k ^\T k \{l + e 2 ^ k ). (51) 

We then obtain 

\A k \ 2 ~ |T fc | 2 cos 2 Xfc (52) 

For potentials characterized by large deviations from parity symmetry, the transmission probability 
\Ak\ 2 may differ significantly from \T k \ 2 . In fact, for Xk — f > the transmission probability vanishes 
even if \Tk\ 2 is non-zero. 



3.4 Delay time 

Next, we examine the case of an initial state localized around and concentrated around 

the value k = p for momentum. In particular, we consider a wave function ipo(k) = Uo(k — p)e lhx °, 
where Uo(k) is a positive- valued, square- integrable function centered around k — 0. Then, the 
integrand in Eq. f HHj) involves a rapidly oscillating phase exp[ik(xo + L) — iE k t + 6k], where 
6k = ImlogAk- Often, but not always, this implies that the integral is strongly suppressed unless 
the phase is stationary for k — p, i.e., unless 

x + L — v p t + 6' p = 0. (53) 

In general, the stationary phase approximation is valid if the probability density Eq. (|48p has a 
single peak. In this case, the solution to Eq. (153"]) . t{p) := (x + L + 8') / v p , determines the location 
of the peak. 

In absence of the barrier, and for the same initial state, the time-of-arrival probability distribu- 
tion is peaked around to(p) = (xq + L)/v p . Thus, when comparing two time-of-arrival experiments, 
one with the barrier and one without, particles in the former are detected with a delay 

Uip) = t-io{p) = d' p /v p = —Im 

when compared with detected particles in the latter. For parity symmetric potentials, Ak = Tk 
and the delay time Eq. (15^1) . when defined, coincides with the Bohm-Wigner time [26]. In general, 
they differ. 

Given Eq. fl54|) one may infer that the time particles spent inside the barrier region equals 

T p = t d (p) - d/v p . (55) 

However, this inference is based on an analogy from classical physics and it does not follow 
from the rules of quantum theory. The time delay t^p), Eq. fl5^|) . is obtained from the comparison 
of probability distributions from different experiments; it is not directly measured in a single 



9 log A k 
dk 



_ k=n 



(54) 
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experiment. The delay time is not a quantum observable or a random variable of the theory. It 
is a temporal parameter that appears in the probability distribution P(L, t) for a specific class 
of initial states. This parameter may incorporate significant information about properties of the 
barrier, but only if it captures a significant physical feature of the probability density P(L,t). 

Moreover, Eq. (|54|) follows from a saddle-point approximation to the probability density 
P(L,t), Eq. (jJSJ) that is valid only if the integrand in Eq. (148 j) has a single maximum. This 
is not true in general, as we will see in the study of the double square barrier in Sec. 4.3. In 
such cases, the naive application of the saddle point approximation leads to erroneous physical 
conclusions. In general, higher moments of the probability distribution P(L, t) contain significant 
information and the temporal aspects of the tunneling process cannot be adequately described in 
terms of a single parameter such as the delay time. 

In what follows, we will refer to r p , Eq. ( 1551) . as "tunneling time", for convenience, without 
committing to its interpretation as the transit time of the particle through the barrier region. We 
will view t p as as a useful time parameter that characterizes the probability distribution P(L,t). 

4 Application: Piecewise constant potentials 

In this section, we construct the probability distribution for the time of detection for two important 
examples of tunneling potentials: the square barrier and the symmetric double square barrier. Both 
potentials are parity-symmetric, so the amplitude A^, Eq. (jHJ), coincides with the transmission 
amplitude T k . 

We consider an initial state i[>o(k), centered around X — — Xq and well concentrated around a 
mean momentum k — p, 

Mk) = u (k-p)e lkxo , (56) 

where Uq is an even, positive- valued wave function centered around k = with width a p . We 
assume that the absorption coefficient a(k) does not vary significantly with momentum, so that 
we can set it equal to unity. Then, Eq. (|15|) takes the form 

(57) 

In what follows, we evaluate Eq. fl57|) for the square barrier and the symmetric double square 
barrier. 

4.1 Junction conditions 

Explicit calculations of tunneling time are usually performed for piecewise constant potentials. 
In such potentials, eigenstates of the Hamiltonian are identified from junction conditions on the 
points of discontinuity. The Hamiltonians hi and /12, Eq. (TTOT) are non-local operators, and their 
junction conditions turn out to be different from the ones usually considered in the literature. 
Here, we present the derivation of the junction conditions. 

We consider the Hamiltonian hi, with a piecewise constant potential. For piecewise constant 
potentials, V(x) = V(x) except for the points of discontinuity. Hence, except for the points of 
discontinuity, the eigenfunctions of hi are of the form Ae tkx + Be~ lkx where k may be either real 
or imaginary. Let us assume that x = is a point of discontinuity, and consider a neighborhood 
U around x = with no other point of discontinuity. Let g(x) be an eigenstate of hi with energy 



P(L,t) 



f°° dk r. 



v k \T k u (k - p)e 



ik(L+xo)—iEkt 
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E. Assuming x e U, g"{x) = —k 2 g(x) for x < 0, and g"(x) = —k%g(x) for x > 0, where /ci and £; 2 
are real or imaginary constants. Then, 

dxhi^(x) = E / dxi/j(x), (58) 



for some constant e that will be eventually taken to zero. For continuous g(x), all terms in the 
above equation that do not involve derivatives vanish at the limit e — > 0. Hence, as e — > 0, 

J dx^/-d 2 + m 2 g(x) + ^ dx<5V#(x)) = 0, (59) 

where 5V = [[V, h l J\ h~ 1/2 ]. 

We write \J p 2 + m? = Yl^=o a nV 2n i m terms of the coefficients a n obtained from the binomial 
expansion. The first term in Eq. (J5"9~]) becomes 

n=l n=l 

- an m 2n - 2 9 \e) - ( fcl )*-V(- e )) = - ^l 2 n 2 ^(e) + ^ " 1 g'(-e) = 

/Co fCi 

n=l z 1 

- ^ p = m 9 '(e) + ^ + f - m g'(-e). (60) 

k 2 k 1 

Writing 1/ ^p 2 + m 2 = b n p 2n , the operator 5V can be expressed as a series 

2n 2m 

^ = ^oE a «E & ™EE^ d '^ 2n+2m ""> ( 61 ) 

n m k=0 1=0 

where Vq is the potential jump at x — 0. Consider the integral of a single term in the series above 
J € dxp k+l 5'(x)p 2n+2m ~ k - l g(x) = (- t ) 2n + 2 ™ f dx + 1 ^ 5 (1+r) (x)9 2n+2m "^(x) = 

= (-i) 2 "+ 2 ™ £ ( k + / ) (-1)- f dx5'(x)d^+ 2m g(x). (62) 

The sum over r involves the term Ylt=o ( ^ ^ J ( — = 0- ~ ^Y = 0- Hence, the integral above 

vanishes, and so does the second term in Eq. ( l59l) . 
Our final result is the junction condition 

Vm 2 + k 2 -m ^m 2 + kl-m , 

U 9 (°+) = u 9 (0-), (63) 

which is to be implemented together with the continuity condition 

9(0+) = g{0-) (64) 
For \ki\ « m and |/c 2 | << m , Eq. f )63|) reduces to the non-relativistic junction condition g'(Q+) = 
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4.2 Square barrier 

4.2.1 Transmission amplitude 

Next, we consider the square barrier potential, 

^) = {^;^[:|| (65) 

For E — Vq < m, we apply the junction conditions (1631) at x = ±~, to obtain the transmission and 
reflection coefficients 

g— ikd 

k cosh(Afc<i) — irjk sinh(Afccf) ^ ^ 

. e" ife Vfc sinh(A fc d) 
cosh^AfcCtj — ir]k sinh(Afca) 



where 



and the functions 



Xk = yfirf -(E- V ) 2 , (68) 



Vk = - ( e k - — ) (69) 



2 V e k 

1 ( 1 

- e fc + - 

2 V e k 



Pk = I ( e k + — ) (70) 



are expressed in terms of the quantity 



A fc y/m 2 + k 2 -m 

ek = T / o fr ( 71 ) 

K m — y m z — A£ 

In the non-relativistic limit, = fc/A^ and Eqs. (1^61— [67]) reduce to the standard textbook 
expressions for the square barrier potential. 

4.2.2 Tunneling time 

Expressing the transmission amplitude (|66|) as = |Tfc|e^ fc , we note that the phase <pk varies much 
faster with k than \T^\. Hence, when computing the probability density Eq. (1571) for a sufficiently 
narrow initial wave-packet, we may assume that \T^\ ~ \T P \. The velocity varies slowly with k 
in comparison to the phases, hence, we also set v k = v p . Expanding the phase factor to first order 
in k — p, 4>k = 4> p + (p' k (k — p), we obtain 

P(L,t) =v p \T p \ 2 \u [v p t-(L + x + (f)' p )}\ 2 , (72) 

where uo(x) = f ^e~ lkx uo(k) is the wave function uq in the position representation. The proba- 
bility density Eq. (1721) exhibits a single peak at t = [L + xo + (fi' p ), thus leading to the identification 
of delay time td{p) = <$>' p , in accordance with Eq. (1541) . Using Eq. ( 1661) . we obtain 

t d (p) = -- + r p , (73) 
v p 
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Figure 1: Tunneling times r p as a function of the incoming momentum p G [0, -y/ 2mVo + V^ 2 ] , for different 
values of the potential Vq and barrier length d = 5000m -1 m, m the particle's mass. 



where the tunneling time is 



r _ ~ ^ m \ P X " d + m ^(p + ^) sinh ( V) cosh(A p rf) 
cosh 2 (A p (i) + 7] 2 sinh 2 (A p c?) 

Eq. applies to all values of p corresponding to tunneling: < p < \/Jrn + Vq) 2 — m 2 . 
For a long barrier e~ Xkd « 1, 

„-ifcd„-A fc d 

T fc = 2- : , (75) 

1 - lVk 



and the tunneling time is independent of d 



Pp A , 1 



r p = m-f4(- + -)• (76) 
l + rjp P K 



4.3 Double square barrier 

Next, we consider the case of the symmetric double square barrier, that is, the piecewise constant 
potential 

VM _ / xe (-oo,-a-f)U(-|,0U(a+|,oo) , . 

W \Vq xe[-a-§,-|]u[§,a + |] ■ 

This potential corresponds to two identical square barriers of width a separated by a distance r. 
4.3.1 The transmission amplitude 

We compute the transmission probability by implementing the junction conditions ( |63i) at x = 
±(a + |) and at x = ±~. We find, 

g— ifc(r+2a) 

(78) 

(cosh(Afcd) — zr/fc sinh(Afca)) 2 e~ 4fer + pi sinh 2 (Afca)e ifcr ' 

where A&, r\\~ and are functions of k that refer to a single barrier and they are defined in Eqs. 

(EBr-raa. 
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The transmission amplitude Eq. (J78"]l can also be expressed as 

rp2 

Tk = 1 _ K>2 p 2ik{r+a) ' ( 79 ) 

where and i?ofc stand for the transmission and reflection amplitudes for the single barrier, as 
given by Eqs. (|66|) and (!67|) respectively. 

Writing T ok = \T ok \e^ k and .Rofc = -i\Rok\e i<i>k , Eq. (JTSJ becomes 



T 2 

0/c (80) 



1 + \R 0k \ 2 e 2i ^ r + a ')+^} 



Eq. (180]) implies that for cos [2(k(r + a) + = —1, |T fc | = 1. It follows that any momenta 
< A;„ < a/VJ) 2 + 2m Vo that are solutions to the equation 

2fc n (r + a) + fcn = (2n + l)vr, n = 0,1,2,... (81) 

correspond to resonances of the double-barrier potential. 

4.3.2 Large inter-barrier separation 

Next, we expand the denominator of the transmission amplitude Eq. (|80|) as a geometric series, 

oo 

T fc = |T 0fc | 2 ^(-l) n | J R fc| 2n e 2i[nfe( - R+a)+(n+1) ' /,fc] . (82) 

n=0 

Substituting into Eq. ( 1571) . we find that the probability density P(L,t) involves the amplitude 

OO /»00 

A(L,t) = ^(-l) n / ^^^(fc-^lTofcH^ofcl^e^' ^"^^ 1 "^^^ (83) 

n=0 ■'O 

We evaluate the integral Eq. ( 1831) using the same approximations as in Sec. 4.2.2. We set 
|T fc = \T p\, \Rok\ = \Rap\, Vk = v p , and we expand the phase factor to first order in k — p. Taking 
the limits of integration to k G (— oo, oo), we obtain 

A{L,t) = ^\T 0p \ 2 e ipL - lEpt (84) 

oo 

^(-^"liJopl^e^^^+^^Uot^o + L) — v p t + 20; + 2n{r + v p r p )}., 



oo 

X 

n=0 



where r p is the tunneling time, Eq. ( 155|) . and «o(x) is the wave function u transformed to the 
position representation. 

The amplitude Eq. ( 18 5 p involves a sum of terms labeled by n, with each term peaked around 
time 

t n = [x + L + <f)' p + 2n(r + v p r p )]/v p . (85) 

The distance between two successive peaks is 

At = 2(r/v p + Tp). (86) 

If the position spread a x > l/(2a p ) of the wave-function uq satisfies a x << v p At, then there is no 
overlap between the terms in the sum Eq. (|85|) . Hence, the amplitude Eq. (185]) is close to zero 
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except for neighborhoods of width a x /v p peaked around the values t n , Eq. (1551) . The first peak in 
the amplitude is at time 

t Q = (L + x Q )/v p + 2(f)' p /v p . (87) 

The next term in the series Eq. (1851) has a peak at to + At. The delay by a factor At can be 
attributed to the fraction of particles transmitted through the first barrier, reflected at the second 
barrier, reflected at the first barrier and exiting the second barrier in its second attempt. Similarly 
the peak at time t n corresponds to particles that were transmitted through the first barrier, they 
were subsequently reflected n times at the second barrier and n times at the first barrier, and 
finally exited after their (n + l)-th attempt on the second barrier. The attempts of crossing the 
barrier are statistically independent: a particle exiting in its (n+ l)-th attempt has been reflected 
2n times. Since the amplitude for a single reflection is Ro p , the contribution of the successful exit 
by the (n + l)-th attempt is suppressed by a factor \R 0p \ 2n . 
The detection probability P(L,t) = |„4(L,t)| 2 is 

oo oo 

P(L,t) = v p \T 0p \ 4 J2Yl \Ro P \ 2{n+m) e iin - m)(3p uo(v p At[n - (t - t ) / At])u (v p At[m - (t - t )/At]),(88) 

n=0 m=0 

where 



(3 P = (2p(r + a) + 2<p p + tt) mod (2tt) (89) 

For a x « v p At, the function u in Eq. ( 155|) approximates well a delta function. Delta 
functions satisfy the property S(x)5(y) = 5 ( gO ~ 2/)- We expect that an analogous equality 
also holds approximately for smeared delta functions, i.e., for families of functions F a (x), such that 
lim a ^ F a (x) = 5(x), when a is close to zero. Hence, we can write uo(x)uo(y) = uo[^(x+y)]uo(x — y) 
in Eq. (158"]) ; the equality is exact for Gaussian wave-functions. We also rearrange the double sum 
in Eq. (!88~]) as J2n=o ^2m=-n> where N = m + n and N = n — m. Then, the time dependence of 
Eq. (1551) is contained in a term 



\Ro P \ 2N u (v p At[-N - (t - t )/ At}). (90) 

N=0 

Thus, the probability distribution P(L,t) is an infinite sum of terms, each sharply peaked 
at times t^ = t Q + \NAt and suppressed by a factor \R 0p \ 2N . If the time-scale of observation 
is much larger than At and |T 0p | 2 << 1, we can approximate the probability distribution by a 
smooth function connecting the peaks of Eq. (|90|) . modulo an overall normalization constant. In 
this approximation, the probability distribution P(L,t) vanishes for t < to. For t > to, P(L,t) is 
proportional to l^ol 2 ^ = e^ 10 ^ 1 -^! 2 ) ^ e -^|To P | 2 ; where n = 2(t - to)/ At. 

Hence, 



p ^ = \cJ^\Vlr (91) 

where 



0\T |2 \T |2 

l p~ At -Ji+T- 1 ] 

Vp P 
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Figure 2: The probability density P(L,t), Eq. ([88]) . for particle detection for the double square barrier 
potential, as a function of the time t — to after first detection. The probability density is a sum of peaks 
separated by distance At. The peaks in the distribution decay exponentially with rate F p , Eq. ([92]h 



is a decay coefficient associated to the barrier, and C is a positive constant. 

Eq. ( 1911) applies for (i) p far from a resonance frequency of the double barrier system and (ii) 
for wave-packets with position spread much smaller than v p At. 

The coefficient C Eq. (I9TT) is determined by the requirement that the time-integrated proba- 
bility is consistent with Eq. (!50|) . i.e., the requirement that 





Tq p 


4 


1 + 


Rop 


4 


+ 2 \Ro p 


2 cos (3 P 



P(L, t)dt ~ |T P | 2 = - ■ - |4 ; p — - . (93) 



This leads to 



^ = { iT/rjVw !>£ • < 94 > 

The physical interpretation of Eq. (I94p . expressed in classical language, is the following. The 
first signal arrives at the detector at a time t which corresponds to a delay time td(p) = 24>' p /v p 
that is twice the delay time of the single barrier. The associated tunneling time r 2fe for the particles 
at first detection is 

rf = t d {p) + = 2t p + -. 95 

v p v p 

The classical interpretation of Eq. (195]) is that the particles cross the two barriers independently 
and at time equal to the tunneling time t p of each barrier. In the inter-barrier region, the particles 
propagate freely with velocity v p . Particles crossing the first barrier, are trapped in the inter- 
barrier region, and escape only after multiple attempts to cross the second barrier. The detection 
probability of those later-escaping particles decays exponentially with a rate T p . 



4.3.3 Wave-packets narrowly concentrated in momentum 

In Sec. 4.3.3, we considered the regime a p v p At » 1 that corresponds to an inter-barrier distance 
r much larger than the spread of the initial wave-packet. In the opposite regime, of small inter- 
barrier distance r, the peaks in the amplitude Eq. (|85[) essentially overlap. Thus, it is a meaningful 
approximation to substitute the infinite sum in Eq. (185!) with an integral over a continuous variable, 
using the Euler-MacLaurin summation formula. 
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We consider a Gaussian u (k), 



-4vr\ 1/4 -A 

//,-,( A- ) ( —J e 4CT ?. (96) 



Then, Eq. ([55]) becomes 

r 2\ 1/4 



/ r /.oo _ 2 i 



,(97) 



where to is given by Eq. ( 1871) . At is given by Eq. (18611 and j3 p by Eq. ( 1891 . 

For t < t Q , the amplitude Eq. ( |9"7|) is strongly suppressed. In the vicinity of t = t , small values 
of n dominate in the integral Eq. (197|) . so we can ignore the contribution of the term e~ n ' T ° p ' . 
Then, 



a 2 X 1/2 

P(L,t) = t; p |T 0p | 4 ( ^ 



7T 



2<TpVp/\t 2 



(9* 



where /x p = 2o - ^ At > anc ^ er ^ c s t an ds for the complementary error function 

2 ro ° 



erfc(x) = -^= / e~ t2 dt. (99) 

7T 



The probability density (1981) is strongly peaked around t = to, hence, also in this regime, the 
delay time for the double barriers is twice the delay time of the single barrier. 

At later times (t — t » l/(a p v p )), we can take the limit of integration in Eq. (197|) to — oo, 
and drop the term outside the integral. Then, 

P(L,t) = Ke- Tp{t - t0 \ (100) 

where 

K = y^iei-rt (101) 
Av p a p 

is a constant. 

Hence, following a transient period of order l/(a p v p ) after the first particle detection, the 
probability P(L,t) decays exponentially with rate T p . This qualitative behavior is verified by a 
numerical evaluation of the probability density Eq. (IHHj) in the regime a p v p At ~ 1 — see, Fig. 3. 

4.3.4 Resonant tunneling 

The results of Sees. 4.3.2 and 4.3.3 rely on the assumption that the particles' momentum p is far 
from the resonances of the double barrier, Eq. ( 1811) . Here, we consider the probability distribution 
P(L,t) for an initial state centered on a momentum p that is close to a resonance frequency ho- 
Near the resonance k = k , the transmission amplitude is approximated by 

\T ok \ 2 e 2i< t >k 

Tk = ~ ~ - 2i{k-k ){r+a+<j>' ) ■ ( 102 ) 



1 - |^ofcol 2 e 







For |T fc | << 1, the modulus square |Tfc| 2 is of the order of |T fc | 4 except for a narrow strip of 
momenta around k where \Tk\ 2 approaches unity. These momenta will dominate in the probability 
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P(t) 

0.35 1 
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Figure 3: The probability density P(L, t), Eq. ([88]) . as a function of the time t, in the regime a p v p At ~ 1. 
The probability exhibits a sharp peak at the time to of first detection, followed by a sequence of shorter 
peaks that decay exponentially with rate T p , given by Eq. ([92]) . 



distribution P(L,t). Hence, we expand the phase factor e 2 ^ fc fc °) (r+a+ ^fco* ) ) keeping only first order 
terms with respect to k — k$. The transition amplitude becomes 



1 fc 

where Tk is the decay coefficient, Eq. (1921) . evaluated at the resonance frequency. 
In this approximation, the squared amplitude 



(103) 



\T k \ 2 = ^ • (104) 



is a Lorentzian of width Tk /(2vk ) centered around k = k . 

We consider a initial wave-packet centered around momentum k = p, such that p lies within 
the peak of the Lorentzian Eq. (I104p . The probability density P(L, t) is obtained from the squared 
modulus of the amplitude 

A(L,t)= J u (k-p)e^^- E ^ (105) 

Jo ^ l-i^-(k-k G ) 

1 fe 

We assume that the momentum spread a p of 0o is much smaller than p, so that we can enlarge 
the integration range to (— oo, oo). We also expand the phase factor to first order in q = k — p, to 
obtain 

/°° tin 1 
— ^ Mk ~ p)e- iyVpit - to) . (106) 
-oo 27Tl_i^JL[ g _ (ko-p)) 

We calculate the integral in Eq. (I106P analytically, by choosing a Lorentzian wave-packet 

u (k) = "^^tV" ( 107 ) 
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Figure 4: The probability density P(L,t), Eq. (|88p . as a function of the time t, near resonance and for 
an initial Lorentzian wave-packet, Eq. (jl07H . The probability exhibits a sharp peak at the time to of first 
detection, and an exponential decay phase with rate Tk , given by Eq. (f92|) . 



Then, we obtain 
A(L,t) = 



ir k Q y^F e ip(xo+L)-iE p t+i<pp 



I _|_ [ fc o~P _ I Ffc o 



— I 



r fc „ 



2f fe o"p 



-«p<Tj,|t— to I 



t < t 







_2i e - iv p( k o-p)( t - t o) a 2 - 1 * — *o I _|_ 



fco— p 



- k 



2v kQ a p 



+ 1 



-v p a p \t-t \^ t>tc 



:io8i 



We note that the amplitude is suppressed for times t < t . The first detection signal appears 
around t = t . Again, this implies that the time delay due to the barrier is 2(j)' p /v p . 

The probability density P(L,t) = \A(L,t)\ 2 is computed straightforwardly. In the regime of 
exact resonance (\ko — p\ « a p ), and for Tk >> a p v p , it simplifies 



P{L,t) 



-v p a p \t—to\ 
k 



t < t 



a p V k [ 4e" rfc o( i - i o) + e -2a p v p (t-t ) + Ae -(^L + a pVp )(t-t ) cos [ Vp (fc - p)(t - *„)], t > t . 



(109) 



The instant to of first detection is followed by a transient regime, and at times t such that t — to » 
(o"pfp) _1 , the probability decays purely exponentially with rate Tk - 

Hence, also in the case of resonance the double square barrier is characterized by two time- 
scales: the delay time 2(p' p /v p and the decay rate T ko of the detection probability at later times. 



4.3.5 Multiple resonances 

For large inter-barrier distances, the difference between two successive resonance momenta is of 
the order of n/r. Hence, for wave packets with significant momentum spread o p ~ n/r, two or 
more resonances may contribute significantly to the detection probability. In this case, we express 
the transmission amplitude as 



T, 



1 



— I 



(/u ky> 



(110) 



where k n correspond to the resonance frequencies. 



22 




Figure 5: The probability density P(L,t), Eq. (|88p . as a function of the time t, for an initial Lorentzian 
wave-packet, Eq. (I107p . that overlaps with two momentum resonances ki and k<i- The probability exhibits 
a sharp peak at the time to of first detection, and an exponential decay phase with rate ~ T^, 
modulated by oscillations of frequency ~ v p \k2 — k\\. 



Using the approximations employed in the derivation of Eq. (1106} we obtain 

A(L,t) = ^e^-^Y: jZt l- {kn - P )) U0{k -^~ iqVp(t ~ t0) - ( U1 ) 

The amplitude A(L, t) is a sum of partial amplitudes A n (t, L), each incorporating the contribution 
from a single resonance. For the Lorentzian initial state q. (I107p . each amplitude A n {L,t) is of 
the form ( 11081) . The resulting probability distribution has the same structure as Eq. ( 11081) around 
t = to, since to depends only on the mean momentum of the wave-packet. However, at later times 
the behavior differs. To a first approximation, the decay constants at different channels are equal 
Tk n — T p . Then, the dominant contribution to P(L,t) at later times is 

P{L,t) ~ | e -^( fe ™-f)(*-*')|2 e -r P (t-t )_ ( 112 ) 

n 

Thus, at later times, the exponential decay of the probability density P(L, t) is modulated by 
oscillations of frequencies v p (k n — k m ) ~ v p /r. 



4.3.6 Summary 

We calculated the probability density P(L, t) for the double-square barrier in different regimes and 
using different approximation schemes. In all regimes, P(L,t) has the same qualitative behavior: 
it is characterized by a sharp peak at the time to of nrs t detection and (sometimes after a transient 
period) then decays exponentially to zero with a rate T p . The probability distribution is determined 
by two parameters that are denned solely in terms of the double barrier's characteristics and the 
incoming particle momentum: the delay time td(p) associated to first detection, and the decay 
rate T p . The delay time td{p) is twice the delay time for a single barrier. From the value of the 
delay time, one is led to a classical picture of tunneling, according to which the particles crosses 
the two barriers independently, each crossing taking time r p each, while propagating freely in the 
inter-barrier region. 

The exponential decay in the probability density P(L,t) at late times is due to the fraction of 
particles that cross the first barrier, and then they are trapped in the inter-barrier region. For these 
particles, the double barrier acts as a potential well, which they can only exit through tunneling. 
The exponential decay law at later times is of the same origin as the exponential decay of unstable 
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states through tunneling. Indeed, its derivation in Sec. 4.3.2 is rather similar to the classic studies 
of a-decay by Gamow [32] and Gurney and Condon [33] — see also, Ref. [20J. 

A key point that emerges from this analysis, is that the temporal aspects of tunneling cannot 
be captured by a single variable. The full construction of the probability distribution P(L,t) is 
necessary. In particular, if one attempts to describe tunneling in terms of the time-delay obtained 
through a stationary phase approximation, as in Eq. (I54p one is led to very different conclusions, 
namely, that there exists a regime, in which the inferred tunneling time is independent of the inter- 
barrier distance r, a property referred to as the "generalized Hartmann effect" [281 129]. However, 
as shown in Ref. [31], a simple stationary phase approximation does not capture the properties 
of the Schrodinger evolved wave function in the double barrier system. A careful implementation 
of the stationary phase approximation in the amplitude A(L, t) of Eq. ( 185]) leads to the results 
presented here. 

5 Conclusions 

In this work, we developed an approach to tunneling time based on the definition of probabilities for 
quantum temporal observables. The tunneling time is defined in terms of a probability distribution 
P(L, t) associated to time-of-arrival experiments in a tunneling set-up. The probability distribution 
P(L,t) contains all information about the temporal properties of the tunneling process. 

We presented a methodology for constructing such probabilities, that applies to any quantum 
system. In particular, we applied this method to spinless relativistic particles described by a com- 
plex scalar quantum field. We derived a simple expression for the probability distribution P(L,t), 
in which all information about the potential barrier is encoded in a complex coefficient con- 
structed from the transmission and reflection coefficients of the potential, and Rk respectively. 
We found that the total transmission probability is proportional to \Ak\ 2 rather than to |Tfc| 2 and 
that this difference may be very pronounced for highly asymmetric potentials. 

We evaluated the probability distribution P(L,t) explicitly for piecewise constant potentials. 
The case of the double barrier potential is of particular significance, because it provides an explicit 
demonstration that a single parameter (like the tunneling time) does not suffice to describe all 
temporal aspects of the tunneling process. 

In the present paper, we specialized to the standard case of a potential barrier that is defined 
in terms of a classical, background electric field. However, the applicability of the method is not 
restricted to this case. The probability distribution P(L, t) is defined in terms of the two-point 
correlation function of any field theory with no restriction on the dynamics. Thus, the method 
presented here can be applied to a broader set of processes, for example, tunneling in open quantum 
systems (including the effect of dissipation and noise), or to incorporate the backreaction of the 
microscopic particles on the potential barrier. 
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